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Spectroscopic  Emission  Measurements 
of  a  Pulsed  Plasma  Thruster  Plume 


Thomas  E.  Markusic  and  Ronald  A.  Spores 

US  Air  Force  Electric  Propulsion  Laboratory.  Edwards  AFB,  CA 


The  optical  emission  spectra  of  plasmas  produced  by  a  Teflon  (polytetrafluoroethylene  (PTFE))  solid 
'  propellant  pulsed  plasma  thruster  (XPPT-1)  were  measured  and  analyzed.  A  broad  temporally  and 
spatially  integrated  survey  of  the  emission  from  3500  to  7500  A  is  reported.  It  is  shown  that  in  the 
,PPT  discharge  energy  range  surveyed  (15  to  45  J),  the  species  formed  do  not  vary,  that  is.  no  new 
species  are  formed  as  the  discharge  energy  is  increased.  Furthermore,  relative  line  ratios  do  not  vary 
significantly  with  energy,  suggesting  that  the  bulk  thermodynamic  properties  of  the  discharge  plasma 
.  may  not  be  affected  by  changing  the  initial  stored  capacitor  energy.  A  spatially  integrated  but  time- 
resolved  measurement  was  carried  out  which  shows  the  single  species  emission  as  a  function  of  PPT 
•discharge  time.  The  magnitude  of  the  emission  was  found  to  track  closely  with  the  PPT  current  and 
depend  on  the"- discharge  energy.  Relative  line  emission  intensities  are  used  to  determine  the  degree 
of  excitation  equilibrium  and  establish  bounds  on  the  plasma  temperature.  Spatially  and  temporally 
resolved  images  are  used  to  estimate  the  plasma  streaming  velocity  from  time-of-fiight  and  Doppler-shift 
techniques. 


I.  INTRODUCTION 

There  is  presently  a  strong,  renewed  inter¬ 
est  in  Pulsed  Plasma  Thrusters  (PPT)  for  a  wide 
range  of  space  missions.  MightySat  II.1[1].,  sched¬ 
uled  for  a  2000  launch,  will  fly  a  PPT  on-board 
that  is  designed  to  perform  an  orbit  raising  ma¬ 
neuver  from  an  initial  Space  Shuttle  altitude  of 
200  to  215  nm  up  to  greater  than  250  nm;  this  ma¬ 
neuver  will  extend  the  life  of  the  satellite  from  two 
months  to  one  year.  The  EOT  satellite,  sched¬ 
uled  for  a  1999  launch,  will  be  using  a  PPT  for 
attitude  control[2].  PPTs  have  the  potential  of 
fulfilling  the  attitude  control  role  on  a  satellite  at 
greatly  reduced  mass  and  cost.  PPTs  are  also  be¬ 
ing  considered  for  constellation  maintenance  for 
such  missions  as  interferometric  imaging  of  the 
Earth  from  space  or  deep  space  from  an  Earth 
orbit.  The  benefits  of  PPTs  are  its  very  small 
impulse  bits  (~  10-6  Ns)  for  satellite  motion,  re¬ 
liability,  high  specific  impulse  and  overall  ease  of 
system  integration.  These  are  extremely  simple 
devices  with  only  one  moving  part,  the  negator 
spring. 

Even  though  there  are  many  missions  where 
PPTs  are  advantageous  with  existing  performance; 
the  range  of  applications  is  expected  to  greatly 
increase  if  their  overall  performance  can  be  im- 
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proved.  There  are.  two  broad  categories  of  ineffi¬ 
ciencies  in  pulsed  plasma  thrusters  that  limit  the 
present  performance.  These  include  energy  in¬ 
efficiencies  which  pertain  to  the  efficient  trans¬ 
fer  of  stored  capacitor  energy  into  acceleration 
of  the  ionized  propellant  and  propellant  ineffi¬ 
ciencies  which  include  particulate  emission  and 
late-time  vaporization[3j.  The  work  described  in 
this  paper  employs  emission  spectroscopy  to  bet¬ 
ter  understand  the  ablation  plasma  formation  and 
subsequent  acceleration  mechanisms  in  present 
PPTs  in  order  to  improve  overall  future  perfor¬ 
mance.  Both  spatial  and  temporal  resolution  are 
used  to  record  the  emission  spectrum  produced 
during  the  arc  discharge  in  XPPT-1,  a  rectangular 
geometry  thruster  described  in  an  earlier  work[3]. 
Investigated  in  this  work  are  the  identification  of 
constituent  species,  temporal  dependence  of  emis¬ 
sion,  electron  temperature  measurements,  and  ex¬ 
haust  velocity.  _ 

While  analysis  of  the  observed  PPT  emission 
is  useful  in  the  determination  of  plasma  param¬ 
eters,  partitioning  of  energy,  and  the  speed  of 
streaming  motion,  it  also  provides  a  starting  point 
for  more  advanced  diagnostics,  such  as  laser-induced 
fluorescence,  by  identifying  electronic  transitions 
that  are  active  in  the  unperturbed  plasma.  In 
the  present  work  we  describe  two  emission  ex- 

Approved  for  public  release;  distribution  unlimited. 


periments  which  are  primarily  qualitative,  and 
two  that  establish  rough  quantitative  bounds  on 
the  plasma  temperature  and  streaming  velocity. 
The  first  experiment  determined  the  atomic  and 
molecular  species  present  at  any  time  during  the 
PPT  discharge.  The  second  experiment  showed 
how  the  production  of  a  species  depends  on  the 
instantaneous  discharge  current.  The  third  ex¬ 
periment  establishes  bounds  on  the  plasma  tem¬ 
perature  by  the  analysis  of  the  distribution  of  up¬ 
per  states  of  a  fixed  specie  and  degree  of  ioniza¬ 
tion.  The  last  experiment  establishes  bounds  on 
the  streaming  velocity  of  the  accelerated  propel¬ 
lant.  This  is  accomplished  by  imaging  a  side  view 
of  the  emission  gated  at  several  times  during  the 
discharge.  A  second  approach  compares  the  si¬ 
multaneously  measured  front  and  side  view  spec- 
-tral  lines  to  determine  a  Doppler  shift  in  the  axial 
emission  from  which  the  axial  streaming  velocity 
follows  directly. 

The  PPT  PTFB  plasma  is  formed  by  the  arc 
ablation  and  subsequent  dissociation  of  the  C2F4 
monomer.  The  species  observed  in  the  present 
■study  were  FVF  ,  C+,  C++,  and  Ca-  Other 
species,  such  as  C,  were  probably  also  present; 
however,  these  species  do  not  have  any  strong 
lines  in  the  wave  number  range  surveyed. 


II.  EXPERIMENT 

In  this  section  the  details  of  the  equipment 
used  in  all  of  the  experiments  is  described;  then,  a 
more  detailed  description  of  the  four  experiments  1 
performed  is  given.  A  top  view  of  the  experimen¬ 
tal  arrangement  is  shown  in  Figure  1.  The  vac¬ 
uum  chamber  has  three  optical  diagnostic  access 
windows.  It  uses  a  single  turbomolecular  pump 
backed  by  a  mechanical  rough  pump  to  achieve  a 
base  pressure  of  about  1  x  10~°  Torr. 

External  optics  were  configured  to  allow  imag¬ 
ing' of  either  the  front  or  side  view  (or  both  simul¬ 
taneously)  of  . the  thruster.  The  lenses  were  <p 50 
mm  with  f/S  which  provided  an  image  magnifi¬ 
cation  of  2.5.  Periscopes  were  used  to  rotate  the 
imaging  plane  and,  in  the  case  of  simultaneous 
imaging  of  the  front  and  side  view,  to  allow  spa¬ 
tial  separation  of  the  images  on  the  CCD  array. 

A  Czerny-Turner  mount  1.26  m  SPEX  1269 
spectrometer  was  used.  Two  different  diffraction 
gratings  were  employed:  a  1200  g/mm  ruled  grat¬ 
ing  blazed  at  1/n  and  a  2400  g/mm  unblazed  holo¬ 
graphic  grating.  With  the  1200  g/mm  grating 
installed,  the  spectrometer  gives  a  first  order  exit 
plane  linear  dispersion  of  about  6.25  A/mm.  Typ¬ 
ical  slit  widths  were  6-10 /i. 

The  spectrometer  output  was  detected  by  ei¬ 
ther  a  Princeton  Instruments  CCD  camera  (model# 
576G/RBE)  or  a  Hamamatsu  photo- multiplier  tube 


VACUUM  CHAMBER 
(-10'5  Torr) 


Figure  1:  Schematic  of  apparatus  for  PPT  emis¬ 
sion  spectroscopy. 

(PMT)(model#  921).  The  camera  has  a  thermo- 
electrically  cooled,  gated,  and  image  intensified 
detector  composed  of  a  two-dimensional  array 
(576x384)  of  22 /r  pixels.  The  overall  spectro¬ 
graph/CCD  system  resolution  was  limited  by  the 
pixel  size  of  the  CCD  array;  the  wavelength  res¬ 
olution  was  found  to  be  0.13  A.  About  80  A  of 
the  optical  spectrum  was  recorded  in  a  single  shot 
with  the  1200  g/mm  grating.  The  minimum  gate 
width  was  100  ns.  The  timing  for  the  gate  was 
coordinated  with  the  PPT  discharge  by  a  Stan¬ 
ford  Research  digital  delay  generator.  The  photo¬ 
multiplier  tube  was  used  to  obtain  continuously 
time  resolved  data.  The  response  time  of  the  de¬ 
tector  used  is  approximately  2.2  ns.  The  spec¬ 
trometer/CCD /PMT  system  spectral  sensitivity 
was  calibrated  using  an  Oriel  quartz-tungsten- 
halogen  (QTH)  lamp  (model#  6334)  and  its  pub¬ 
lished  spectral  irradiance[4]. 
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A.  SURVEY  OF  SPECTRAL  LINES 

An  initial  experiment  was  conducted  to  de¬ 
termine  the  species  present  in  the  plume  at  any¬ 
time  during  the  discharge.  The  spectrograph  was 
scanned  in  80  A  increments  (the  1200  g/mm  grat¬ 
ing  was  used)  from  3500  to  7500  A.  The  optics 
were  arranged  in  the  front  view  path  orientation 
illustrated  in  Figure  1.  A  central  horizontal  slice 
(approximately  2.5  cm  wide)  of  the  PTFE  propel¬ 
lant  face  was  imaged  onto  the  spectrograph  en¬ 
trance  slit  (7 fj,  width)  and  subsequently  onto  the 
CCD  array.  Since  the  XPPT-1  electrodes  are  2.5 
cm  wide,  this  imaging  assured  that  the  discharge 
arc  would  be  captured  by  the  camera  regardless  of 
.  where  it  initiated  or  terminated.  Typical  XPPT- 
1  discharges  range  from  fifteen  to  twenty-five  mi¬ 
croseconds  in  duration,  depending  on  the  discharge 
energy.  The  CCD  gating  was  adjusted  such  that 
the  CCD  shutter  was  opened  approximately  5  /is 
before  the  discharge,  and  closed  approximately  75 
ps  after  the  discharge  terminated  (i.e.  there  was 
no  current  flowing).  Thus,  the  acquired  images 
represent  a  Spatially  and  temporally  integrated 
accounting  of  species  present  during  and  imme¬ 
diately  after  the  discharge.  Scans  were  performed 
at  three  different  PPT  discharge  energies:  15  J, 
25  J,  and  45  J. 

B.  TIME  RESOLVED  SPECTRA 

Once  the  strong  lines  had  been  identified,  time 
resolved  measurements  were  conducted.  The  same 
front  view  optical  arrangement  as  in  Experiment 
A  was  used:  however,  the  CCD  camera  was  re¬ 
placed  with  an  exit  slit  and  a  photo-multiplier 
tube.  For  each  measurement,  the  monochroma¬ 
tor  wavelength  was  set  to  the  peak  intensity  of 
a  previously  identified  transition.  The  entrance 
and  exit  slits  were  set  to  8 p  and  6/j,.  respectively. 
The  PMT  and  PPT  currents  were  measured  si¬ 
multaneously  and  displayed  and  recorded  on  a 
Tektronix  digital  oscilloscope.  Both  the  PMT  and 
the  PPT  were  triggered  by  a  Stanford  Research 
digital  delay  generator.  This  procedure  was  car¬ 
ried  out  for  several  lines  of  each  specie.  These 
measurements  were  also  performed  at  three  dif- 
ferent-PPT  discharge  energies:  15  J.  25  .J,  and  45 


C.  PLASMA  TEMPERATURE 

The  correlation  of  relative  line  intensities  of 
members  of  the  same  species  with  plasma  tem¬ 
perature  was  used  to  determine  an  estimate  of  the 
electron  temperature.  Only  the  carbon  ion  proved 
to  have  a  sufficient  distribution  of  lines  to  be  used 
for  such  an  analysis;  seven  C+  lines  were  obtained 
using  the  CCD  camera.  The  CCD  shutter  gate- ' 
open  time  was  set  to  100  ns  and  the  gate  delay 


was  adjusted  to  make  the  CCD  shutter-open-time 
correspond  to  peak  PPT  current:  the  two  wave¬ 
forms  were  monitored  on  the  digital  oscilloscope. 
The  side-view  optical  setup  was  chosen  to  avoid 
non-uniform  line  distortion  from  the  Doppler  ef¬ 
fect.  The  1200  g/mm  diffraction  grating  was  used 
and  the  entire  optical  path  was  calibrated  for  in¬ 
tensity  using  the  the  QTH  lamp. 

D.  PLASMA  STREAMING  VELOCITY 

Two  approaches  were  used  to  estimate  the  speed 
of  the  plasma  leaving  the  thruster  during  the  arc 
discharge  phase. 

The  first  technique  used  time-of-fiight  to  track 
the  position  of  the  plasma  front  as  it  expanded 
from  the  solid  PTFE  face.  The  emission  was  im¬ 
aged  in  the  side  view  orientation  (see  Figure  1). 
The  CCD  camera  shutter  gate-open-time  was  100 
ns;  the  gate  delay  was  adjusted  to  give  “snap¬ 
shots"  of  the  plasma  at  various  times  during  the 
discharge  cycle.  The  resulting  images  then  show 
plasma  front  position  as  a  function  of  time,  that 
is,  plasma  speed  transverse  to  the  viewing  direc¬ 
tion.  The  CCD  camera  was  rotated  ninety  degrees 
to  allow  the  maximum  amount  of  space  to  be  im¬ 
aged  (i.e.  the  578  pixel  direction  was  aligned  ver¬ 
tically  with  the  spectrograph  entrance  slit(25  //•)).' 
The  spatial  dimension  of  the  CCD  array  was  cal¬ 
ibrated  using  a  Mercury  pen  ray  lamp  masked  to 
a  known  dimension.  The  lamp  was  placed  in  the 
thruster  position  and  imaged  back  to  the  spectro¬ 
graph  entrance  slit  and  then  acquired  by  the  CCD 
camera  allowing  the  number  of  CCD  pixels  per 
unit  distance  to  be  calculated.  It  was  found  that 
approximately  seven  centimeters  of  axial  distance 
in  front  of  the  solid  PTFE  face  could  be  imaged. 
Representative  lines  of  F+,  C+,  and  C++  were 
recorded  at  various  discharge  times  for  later  anal¬ 
ysis.  All  of  the  images  were  acquired  with  25  J 
PPT  discharge  energies. 

The  second  approach  resolved  the  Doppler  shift 
of  the  emitted  radiation  from  particles  moving  ax¬ 
ially  away  from  the  thruster.  Both  the  front  and 
side  views  of  the  thruster  were  imaged  simulta¬ 
neously.  This  obviated  the  need  for  exact  wave¬ 
length  calibration  since  the  side-view  line  was  as¬ 
sumed  to  have  zero  Doppler  shift,  and  therefore 
provided  an  intrinsic  reference  point  for  the  shifted 
front-view  line.  An  Inconel  beam  splitter  was 
used  to  image  both  views  onto  the  entrance  slit 
(11  p).  The  2400  g/mm  grating  was  implemented 
to  provide  maximum  dispersion,  and  the  CCD 
camera  was  in  the  standard  orientation  (longer  di¬ 
mension  in  the  wavelength  direction).  The  CCD 
gate-open  time  was  100  ns, -and  the  images  were 
acquired  at  peak  current  (25  J  discharge  energy). 
Due  to  the  narrow  slit,  short  gate  time,  and  losses 
in  the  beam  splitter,  only  the  carbon  ion  had  suf¬ 
ficient  intensity  to  be  resolved. 
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III.  RESULTS 


A.  SURVEY  OF  SPECTRAL  LINES 

Five  species  were  identified  by  comparing  the 
acquired  emission  lines  with  a  standard  reference^] 
F,  F+,  C+,  C++.  and  C2.  Figure  2  shows  these 
lines;  prominent  lines  are  labeled  to  identify  the 
corresponding  specie.  Some  of  the  lines  are  over¬ 
lapping  multiplets  and  therefore  represent  aver¬ 
ages  in  intensity  and  wavelength.  The  figure  con¬ 
tains  plots  for  three  different  discharge  energy  lev¬ 
els:  15  J,  25  J,  and  45  J  (the  species  are  la¬ 
beled  only  in  the  15  J  case).  All  of  the  plots 
have  been  calibrated  to  remove  instrumental  ef¬ 
fects  that  change  with  wavelength  or  environmen¬ 
tal  conditions:  intensity,  wavelength,  and  back¬ 
ground  signal  (primarily  thermal) .  The  plots  rep¬ 
resent  averages  over  three  discharge  cycles. 

Table  1  gives  the  wavelength  and  term  desig¬ 
nations  (as  found  in  Moore[6j)  of  the  upper  and 
lower  electronic  states  of  the  dominant  transitions 
(i.e.  those  with  the  greatest  emissivities).  Also 
apparent  in  the  plots  are  the  rotational  and  vibra¬ 
tional  bands  of  molecular  carbon— the  so-called 
Swan  bands  [7]. 

The  three  plots  in  Figure  2  are  structurally 
very  similar.  The  only  apparent  difference  is  a 
modest  increase  in  relative  fluorine  emission  at 
higher  PPT  discharge  energies.  Since  a  change  in 
relative  emissivity  is  directly  related  to  a  change 
in  relative  species  population,  the  figure  shows 
that  increasing  discharge  energy  leads  to  a  slight 
increase  in  the  relative  population  of  F+  while  the 
C+  and  C++  lines  show  little  variation  (note:  the 
ionization  potentials  of  C+,  F,  and  C  are  24.38 
eV.  17.42  eV.  and  11.26  eV,  respectively^]).  Two 
conclusions  can  immediately  be  reached.  First, 
increasing  PPT  discharge  energy  does  not  lead 
to  (at  least  in  the  energy  range  of  the  present 
survey)  a  significant  production  of  more  highly 
ionized  species  such  as  F++  and  C+++  (the  ion¬ 
ization  potentials  of  F+  and  C++  are  34.97  eV 
and  47.89  eV.  respectively [5]).  Since  the  addi¬ 
tional  discharge  energy  is  not  going  into  the  pro¬ 
duction  of  more  highly  ionized  species,  it  must  be 
accounted  for  elsewhere.  One  possiblility  is  that 
increased  stored  energy  leads  to  the  the  ablation 
and  ionization  of  more  material,  rather  than  the 
deposition  of  the  additional  energy  into  a  constant 
amount  of  matter.  A  second  possibility  is  that 
the  degree  of  ionization  (i.e.  the  total  percentage 
of  liberated  PTFE  that  is  ionized)  is  increased. 
These  conclusions  are  further  supported  by  the 
fact  that  the  relative  line  intensity  ratios  between 
members  of  the  same  specie  do  not  change  signif¬ 
icantly  as  the  discharge  energy  is  increased.  This 
implies  that  the  electron  temperature  remains  re¬ 
mains  fairly  constant  regardless  of  the  discharge 


energy. 

Two  species  are  conspicuously  absent  in  the 
spectrum:  C  and  CF.  These  species  are  undoubt¬ 
edly  present;  however,  their  strongest  lines  occur 
in  the  UV  and  IR— regions  that  were  outside  of 
the  present  survey. 

B.  TIME  RESOLVED  SPECTRA 

The  time  dependence  of  emission  was  recorded 
for  C+,  C++.  F,  and  F+  at  three  different  dis¬ 
charge  energies  and  is  shown  in  Figure  3.  The 
specific  lines  used  to  produce  the  plots  were  as 
follows:  4267.26  A  for  C+,  4647.42  A  for  C++, 
6856.03  A  for  F,  and  3847.09  A  for  F+.  These 
lines  were  chosen  because  they  gave  the  best  sig¬ 
nal  to  noise  ratio.  Other  lines  gave  qualitatively 
similar  results  with  reduced  signal. 

Table  1:  Term  and  level  designations  of  dominant 
transitions  in  the  PPT  discharge. 


Specie 

Emission 
Line  [A] 

Lower 

Level 

Upper 

Level 

C+ 

3918.98 

3p  2P?/, 

4s  2S1/2 

C+ 

3920.69 

3p  2P.v2 

4s  2S1/0 

C+ 

4267.00 

3d  2D3/2 

« 2nn. 

C+ 

4267.26 

3d  2D5/9 

4f  2F?/9 

C+ 

5889.77 

3d  2D5/2 

4p  2P°/2 

C+ 

5891.59 

3d  2D3/2 

4p2P3°/2 

c+. 

6578.05 

3s  2S1/2 

3P  ’f  V2 

C4 

6582.88 

3s  2S1/9, 

3p  2P?/2 

C+ 

7231.32 

3p  2P 

3d  2D3/, 

c+ 

7236.42 

3p  2P°/2 

3d  2D5/2 

C++ 

4647.42 

3s  3S1 

3P  3p,: 

C++ 

4650.25 

3s  3SX 

3p  3P? 

F  ■ 

6239.65 

3s  4P3 

3p  JS? 

F  ' 

6348.51 

3s  4P, 

3p4S? 

F 

6856.03 

3s  4P3° 

■3p  4D, 

F 

7037:47 

3s  2P., 

3p  2P2° 

F 

.  7127.89 

3s  2P 2  ' 

3p  2P° 

F 

7754.70 

■3s  2P, 

•3p  2D? 

F4 

3501.45 

3p5Pi 

3d  5Dg 

F4 

3847.09 

3s  5Sf 

3p  5P3 

F+ 

4024.73 

■3s  3S? 

3p  3P-> 

F4 

4103.51 

3p  3P3 

3d  3D3 

F4 

4246.23 

4f  5F3 

3d  5D2 
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Figure  2:  Survey  of  PPT  optical  emission  spectrum  from  3500  to  7500  A  for  three  different  discharge 
energies:  a)  15  J.  b)25  J.  c)  45  J. 
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The  illustrated  data  represent  averages  taken  over 
ten  discharge  cycles.  All  data  was  normalized, 
and  therefore  is  in  arbitrary  units  (i.e.  no  corre¬ 
lation  should  be  made  between  absolute  intensi¬ 
ties  of  different  species).  The  peak  current  in  each 
case  was:  9.75  kA  at  15  3.  12.5  kA  at  25  J,  and 
18.5  kA  at  45  J  discharge  energy. 

All  of  the  data  show  that  the  emission  tracks 
closely  with  the  current.  The  peak  emission  oc¬ 
curs  near  peak  current— lagging  the  current  by 
less  than  a  microsecond.  The  delay  in  peak  emis¬ 
sion  may  be  explained  if  we  consider  the  PPT  as 
a  one-sided  ablation  controlled  arc[8].  The  ob¬ 
served  delay  could  then  be  attributed  to  the  pro¬ 
pellant  transit  time  between  being  liberated  from 
the  solid  PTFE  face  and  being  entrained  in  the 
current  sheet  (i.e.  the  solid  face  and  arc  are  sep¬ 
arated  by  a  thin  vapor  layer). 

From  Figure  3  it  is  clear  that  the  emission  pro¬ 
files  are  affected  by  the  discharge  energy  in  two 
ways.  First,  increasing  the  discharge  energy  re¬ 
sults  in  about  one  half-cycle  of  additional  current 
swing  per  ten  Joules  of  additional  energy  (e.g.  the 
15  J  discharge  terminates  after  one  cycle  while  the 
45  J  discharge  terminates  after  two  and  one  half¬ 
cycles).  The  additional  current  oscillations  pro¬ 
vide  an  opportunity  for  the  plasina  to  “re-ignite” 
and  emit  radiation.  Also,  higher  discharge  ener¬ 
gies  give  higher  peak  currents  and,  in  turn,  higher 
populations  of  emitters. 

General  trends  are  evident  in  Figure  3.  As  the 
current  “rings  down”  the  discharge  supports  de¬ 
creasing  populations  of  the  more  highly,  ionized 
species.  This  is  most  apparent  in  Figure  3d), 
where  the  C+'i‘  emission  is  almost  nonexistent  in 
■  the  second  half-cycle.  However,  C+  and  F  seem  to 
defy  this  trend,  at  least  in  the  first  discharge  cy¬ 
cle.  For  example,  in  Figure  3e)  the  C+  emission  is 
actually  twenty-five  percent  greater  in  the  second 
half-cycle  than  in  the  first  half-cycle,  even  though 
the  current  is  only  sixty  percent  of  its  peak  value. 
This  suggests  that,  in  the  second  half-cycle,  the 
current  carrying  electrons  have  insufficient  energy 
to  produce  a  sizeable  population  of  the  species 
with  higher  ionization  potentials  (F+  and  C4"4"); 
therefore,  a  greater  proportion  of  the  discharge 
energy  is  available  to  create  the  more  weakly  ex¬ 
cited  species  F  and  C+;  thus,  the  emission  of  these 
species  is  maintained  near  peak  levels  in  both  the 
first  and  second  half-cycles. 


C.  PLASMA  TEMPERATURE 


Figure  3:  PPT  emission  as  a  function  of  discharge 
time  and  energy. 


When  the  conditions  of  a  plasma  are  such  that 
local  thermodynamic  equilibrium  (LTE)  applies, 
the  populations  of  the  bound  states  follow  a  Boltz¬ 
mann  distribution^].  It  is  uncertain  whether  the 
the  prevailing  conditions  in  a  transient  discharge 
such  as  the  PPT  reach  a  state  of  thermal  equilib¬ 
rium,  or  if  measured  excitation  equilibrium  neces- 
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sarily  implies  complete  local  thermodynamic  equi¬ 
librium  (CLTE).  In  the  present  work,  the  assump¬ 
tion  of  LTE  is  made  a  priori ,  and  explored  by 
referring  to  the  measured  data. 

Relative  emissivities  of  singly  ionized  carbon 
will  be  used  to  calculate  electron  temperature. 
The  lines  must  be  well  resolved  for  accurate  emis¬ 
sivities  (emn),  and  the  transition  probabilities  (Amn) 
must  be  known.  Since  the  populations  of  the  ex¬ 
cited  states  are  given  by  the  Boltzmann  distribu¬ 
tion,  Equation  1  describes  their  relative  emissiv- 
ity: 


where,  for  a  transition  from  upper  state  m  to 
lower  state  n,  vmn  is  the  frequency,  Em  and  gm 
are  the  energy  and  degeneracy  of  the  upper  state, 
respectively,  k  is  the  Boltzmann  constant,  and  T 
is  the  temperature.  A  plot  of  the  quantity  on  the 
left-hand-side  of  Equation  1  versus  Em  has  a  slope 
of  -1/kT.  Therefore,  the  temperature  can  be  ob¬ 
tained  without  N,  the  total  number  density,  or  2, 
the  partition  function. 

For  experimental  data,  the  application  of  this 
technique  is  only  accurate  when  the  distribution 
of  measured  upper  state  energies  (i.e.  the  spe- 
tral  lines  which  represent  transitions  from  these 
states)  is  comparable  to  the  ambient  electron  tem¬ 
perature.  Of  the  dominant  transitions  listed  in 
Table  1.  only  C+  exhibits  a  sufficient  distribution 
of  upper  states,  and  has  satisfactorily  resolved 
lines,  to  be  useful  in  the  present  analysis.  The  dis¬ 
tribution  of  upper  states  listed  in  Table  2  spans 
approximately  4.6  eV  in  energy.  We  can  there¬ 
fore  legitamately  expect  to  have  the  resolution  to 
analyse  temperatures  ranging  from  1  to  10  eV. 

Seven  well-resolved  transitions  of  C"5"  were  avail¬ 
able.  The  upper  level  degeneracies  and  transition 
probabilities  are  given  in  Table  2.  Only  relative 
emissivities  are  required  and  these  were  obtained 
by  numerical  integration  of  the  CCD  output.  An 
eighth  line  (C+  4267.26)  that  was  not  “well  re¬ 
solved”  is  included.  Since  this  line  dominates 
the  spectrum  (see  Figure  2).  its  inclusion  in  the 
present  analysis  is  necessary.  The  line  is  not  well 
resolved  because  it  overlaps  the  C+  4267.0  line 
with  the  available  instrumental  resolution.  The 
C+  4267.26  line  may  be  included  in  the  tempera¬ 
ture  analysis  by  assuming  that  it  is  in  excitation 
equilibrium  with  the  C+  4267.0.  In  this  case,  the 
relative  emissivities  of  the  two  lines  can  be  found 
from  Equation  1: 

O267.0  _  94267.0  A42S7.0  ^4267.0  ~  ^ 

0267.26  <74267.26  A4267.26  </4267.26 

It  should  be  noted  that  this  separating  of  the  two 
overlapping  carbon  ion  lines  in  no  way  invalidates 


Table  2:  Parameters  for  electron  temperature  de¬ 
termination. 


Specie 

Emission 
Line  [A] 

Em 

(cm-1) 

§ro 

A-mn 

(X108  s-1) 

C+ 

3918.98 

157  234 

2 

0.636 

C+ 

3920.69 

157  234 

2 

1.27 

C+ 

4267.26 

168  978 

8 

2.38 

C+ 

5889.77 

162  525 

4 

0.315 

C+ 

6578.05 

131  736 

4 

0-363 

C+ 

6582.88 

131  724 

2 

0.362 

C+ 

7231.32 

145  549 

4 

0.352 

C+ 

7236.42 

145  551 

6 

0.422 

their  inclusion  in  the  Boltzmann  plot  since  both 
lines  have  the  same  upper  state. 

A  plot  of  the  left-hand-side  of  Eq.  1  using  the 
experimentally  acquired  lines  given  in  Table  2  is 
shown  in  Figure  4,  where  each  of  the  points  rep¬ 
resents  data  averaged  over  five  shots.  A  least 
squares  linear  fit  is  also  illustrated.  The  linear  fit 
indicates  an  electron  temperature  of  1.4A0.2  eV 
(~  .16000  K).  The  experimental  data  does  not  fall 
perfectly  on  the  straight  line,  nor  does  it  deviate 
in  any  ordered  manner;  therefore,  no  strong  argu¬ 
ment  for  or  against  LTE  can  be  made.  However, 
since  the  deviation  from  linearity  is  not  large,  we 
can  safely  assume  that  the  calculated  tempera- 


Figure  4:  Relative  emissivities  corrected  for  de¬ 
generacy  and  radiative  rate  vs.  energy  of  upper 
state. 
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ture  is  representative  of  the  mean  thermal  energy 
of  the  electrons. 

D.  PLASMA  STREAMING  VELOCITY 

The  velocity  of  the  ejected  plasma  varies  with 
position  and  time  in  the  PPT.  In  this  section  we 
determine  the  speed  of  the  propellant  immedi¬ 
ately  after  it  is  liberated  from  the  solid  PTFE 
surface  during  the  first  half-cycle  of  current  flow. 

The  first  approach  used  was  time-of-flight.  Fig¬ 
ure  .5  shows  a  typical  current  trace.  Four  times 
are  indicated  by  the  arrows  on  the  horizontal  axis 
of  the  graph:  times  A,  B,  C,  and  D'(0.1  pis,  1.6 
/xs,  3.68  fj, s,  and  4.68  pis,  respectively)  .  These 
times  correspond  to  CCD  “snapshots”  that  are 
illustrated  in  Figure  6.  The  horizontal  axis  in 
Figure  6  is  the  axial  distance  from  the  propellant 
face.  The  physical  position  and  extent  of  the  PPT 
electrodes  is  is  illustrated  above  the  curves. 

Plots  A.  B.  C,  and  D  show  the  emission  from 
C*  as  the  plasma  expands  downstream  early  in 
the  first  half-cycle  of  the  discharge.  To  determine 
the  rate  of  expansion  (i.e.  the  speed)  we  must  es¬ 
tablish  where  the  plasma  “front”  is  in  each  case. 
This  position  is  somewhat  subjective;  therefore, 
any  reasonable  measure  that  is  consistently .  ap¬ 
plied  to  all  cases  will  be  valid.  We  choose  to  de¬ 
fine  the  plasma  front  as  that  position  where  the 
degrading  emission  reaches  1/e  of  its  peak  value. 
These  points  are  indicated  in  Figure  6  by  the  ver¬ 
tical  dashed  lines  that  intersect  each  curve.  Com¬ 
bining  the  temporal  information  from  Figure  5 
with  the  spatial  information  from  Figure  6  we 
can  compute  the  speed  of  the  expanding  front. 


Figure  5:  Typical  PPT  current  trace  with  CCD 
image  timing  locations,  on  horizontal  axis. 


Comparing  the  emission  at  points  A  and  B  we 
find  v front  8.1  km/s,  B  and  C  give  VfTOnt  ~ 
15.5  km/s.  and  C  and  D  give  Vfront  ~  14.7  km/s. 
Apparently  the  carbon  ions  continue  to  acceler¬ 
ate  between  times  A,  B,  and  C  and  have  reached 
a  terminal  speed  of  approximately  15  km/s  be¬ 
tween  times  C  and  D.  Similar  measurements  were 
made  on  C4"1"  and  F4";  the  speeds  were  found  to 
be  approximately  14.1  km/s  and  15.9  km/s,  re¬ 
spectively. 

A  second  diagnostic  to  determine  streaming 
velocity  is  the  Doppler  shift  of  emitted  radiation. 
The  speed  v  of  an  emitting  particle  moving  rela¬ 
tive  to  a  stationary  observer  is  given  by, 


where  c  is  the  speed  of  light,  A0  is  the  rest  wave¬ 
length.  and  A  is  the  measured  wavelength. 

Both  the  side  and  front  view  of  the  thruster 
were  simultaneously  imaged  onto  the  entrance  slit. 
This  allows  for  an  intrinsic  reference  point  for 
Doppler  shift  measurements.  More  specifically, 
the  measurement  of  the  Doppler  shift  does  not  de¬ 
pend  on  wavelength  calibration  of  the  spectrom¬ 
eter  since  the  side  view  of  the  thruster  provides 
a  simultaneous  (and  presumably  unshifted)  refer¬ 
ence  line.  I 

Figure  7  shows 'the  C+  6578.05  A  line  acquired 
in  the  simultaneous  front/side  arrangement.  Two 
peaks  separated  by  approximately  0.21  A  are  ev- 


Figure  6:  Axial  spatial  distribution  of  emission 
at  four  different  times.  The  PPT  position  and 
electrode  extent  is  illustrated  above  the  graph. 
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Figure  7:  C:  emission  line  representing  simulta¬ 
neous  front  and  side-view  images. 

ident.  From  Equation  3  we  compute  the  parti¬ 
cle  speed  to  be  9.6  ±5.5  km/s,  where  the  margin 
of  error  is  a  result  of  the  maximum  spectrometer 
resolution.  Only  the  carbon  ion  line  had  sufficient 
intensity  to  be  resolved  with  the  narrow  entrance 
slit  and  short  gate  width  used  in  this"  experiment. 

The  Doppler  shift  illustrated  in  Figure  7  was 
measured  at  2.1  y s  (see  Figure  6).  The  time- 
of-flight  estimated  speed  in  this  region  (between 
points  B  and  C)  was  found  to  be  between  8.1  km/s 
and  15.1  km/s.  Therefore,  the  time-of-flight  and 
Doppler  shift  techniques  have  been  shown  to  give, 
within  the  intrinsic  system  uncertainty,  compara¬ 
ble  estimates  of  the  plasma  streaming  speed. 

IV.  DISCUSSION. 

The  primary  goal  of  this  study  was  to  deter-' 
mine  the  constituent  particles  that  are  present  in 
a  PTFE  PPT  discharge,  and  when  these  parti¬ 
cles  form  relative  to  the  current  pulse.  Secon¬ 
darily,  approximate  quantitative  estimates  of  the 
plasma  temperature  and  streaming  velocity  were 
desired.  _  In  this  section,  sources  of  error  m  the 
acquisition  and  interpretation  of  the  experimen¬ 
tal  data  are  explored.  These  sources  of  error  are 
especially  pertinent  to  the  validity  of  the  quan¬ 
titative  sections  of  the  paper..  Without  careful 
application  of  tomographic  techniques  to  resolve 
phenomena  spatially,  emission  spectroscopy  pro¬ 
vides  only  qualitative  information  or,  at  be'st,  al¬ 
lows  one  to  establish  bounds  on  the  value  of  a 
spatially  integrated  experimental  parameter. 

Calculations  were  carried  out  to  estimate  the 
optical  thickness  of  the  plasma.  The  absorption 
coefficient  due  to  inverse  Bremsstrahlung  Kv  is 
given  by  Eq.  4,  where  Z  is  the  average  charge 


(one)  and  v  is  the  frequency  of  the  light  [11] 

Kv  =  3.69  x  10s  (ZNe2)/{T1/2is3)  .  (4) 

Losses  due  to  stimulated  emission  are  assumed  to 
be  insignificant  under  our  conditions  and  the  cor¬ 
rection  term  for  it  has  not  been  included.  At  Ne  = 
1  x  1017  cm-3,  Te  =  16000k,  and  A  =  5500  A,  Eq. 
4  yields  K„  =  1.8  x  10-4  cm-1.  Since  the  plasma 
extent  is  a  few  centimeters,  the  weak  absorption 
(l/Ku  =  55.6  m)  due  to  inverse  Bremsstrahlung 
does  not  affect  our  quantitative  spectroscopy. 

The  data  presented  in  the  spectral  survey  and 
time  resolved  emission  sections  of  the  paper  are 
essentially  qualitative.  The  data  were  compiled 
as  averages  over  several  discharges  and  therefore 
remove  any  non-periodic  phenomena. 

In  reference  to  the  plasma  temperature  section 
of  the  paper,  a  very  important  question  must  be 
posed:  “Where  exactly  in  the  plasma  is  the  elec¬ 
tron  temperature  1.4  eV?”  The  acquired  CCD  im¬ 
ages  represent  spatially  integrated  emission;  the 
collected  light  is  coming  from  many  different  lay¬ 
ers  within  the  plasma  which  undoubtedly  have 
different  temperatures.  "We  make  the  assumption 
that  the  preponderance  of  the  collected  light  is 
generated  from  within  or  very  near  the  arc  col¬ 
umn.  and  the  quoted  temperature  therefore  ap¬ 
plies  only  to  that  region  of  the  plasma. 

The  time-of-flight  data  and  analysis  should  be 
an  accurate  representation  of  the  initial  plasma 
expansion.  The  C+  lifetimes  are  short  in  compar¬ 
ison  to  the  CCD  camera  gate  width.  The  plasma 
translation  is  less  than  0.3  mm  while  the  cam¬ 
era  shutter  is  open.  Two  other  possible  sources 
of  error  exist  in  the  time-of-flight  data.  We  as¬ 
sumed  that  the  regions  that  were  emitting  light 
represented  the  spatial  extent  of  the  gas  that  had 
been  accelerated  from  the  thruster:  however,  if 
processes  such  as  recombination  were  to  deplete 
the  leading  edge  of  the  expanding  front  of  emit¬ 
ters,  we  would  not  observe  emission  from  that 
region.  Consequently,  our  estimate  of  the  speed 
would  be  lower  than  the  actual  exhaust  speed  of 
the  thruster.  Second,  the  emission  from  the  arc 
region  of  the  plasma  may  overwhelm  the  signal 
coming  from  other  regions  of  the  plasma.  This 
would  have  the  effect  of  not  allowing  the  accurate 
determination  of  axial  extent  of  the  plasma.  This 
effect  would  be  especially  important  in  the  early 
stages  of  the  discharge,  such  as  those  indicated 
as  points  A  and  B  in  Figure  6,  and  would  lead 
to  the  calculation  of  the  arc  propagation  speed 
rather  than  the  plasma  streaming  speed.  This 
may,  in  part,  explain  the  difference  in  the  cal¬ 
culated  speeds  between  point  A.B  and  B.C.  The 
Doppler  shift  data  is  not,  by  itself,  an  accurate 
measurement.  Clearly  the  two  lines  are  not  com¬ 
pletely  resolved  (Rayleigh’s  criterion  is  not  satis¬ 
fied).  However,  accounting  for  the  spectroscopic 
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system  resolution,  the  Doppler  measurement  does 
allow  us  to  establish  an  upper  bound  on  the  car¬ 
bon  ion  speed  of  about  15  km/s.  In  short,  while 
the  measurement  does  not  provide  an  accurate 
quantitative  value  of  the  speed,  it  does  establish 
bounds  on  the  speed  which  support  the  time-of- 
flight  findings.  Also,  as  was  mentioned  in  the  tem¬ 
perature  discussion,  the  lack  of  spatial  resolution 
in  the  measurement  implies  that  we  are  measuring 
the  Doppler  shift  of  a  distribution  of  velocities;  we 
are  integrating  emission  from  particles  that  are 
at  the  face  of  the  PTFE  (with  presumably  low 
speeds)  through  particles  that  have  reached  their 
terminal  speeds.  Consequently,  the  Doppler  shift 
data  may  be  dominated  by  emission  from  slow 
moving  particles  near  the  PTFE  surface,  which 
does  not  accurately  represent  the  final  distribu¬ 
tion  of  velocities. 

Line-width  and  line-shape  analysis  is  a  com¬ 
mon  technique  used  to  extract  plasma  parame¬ 
ters.  This  type  of  analysis  was  not  conducted  in 
the  present  study  for  two  reasons.  First,  for  the 
estimated  plasma  conditions,  the  contributions  to 
line  broadening  due  to  Doppler  broadening  and 
the  Stark  effect  are  found  to  be  of  comparable 
magnitude.  Therefore,  it  would  be  difficult  to 
deconvolve  the  two  effects  in  a  given  line.  Sec¬ 
ond,  the  estimated  Doppler  and  Stark  widths  are 
smaller  than  the  available  resolution  of  the  spec¬ 
troscopic  system. 

Measurements  similar  to  those  in  the  present 
work  were  carried  out  by  earlier  researchers[12]. 
In  particular,  exhaust  velocities  were  determined 
using  time-of-flight  and  Doppler  shift  techniques. 
The  range  of  axial  speeds  found  in  that  study 
were:  11.0  to  28.0  km/s  for  C4".  13.0  to  42.0  km/s 
for  C4-4",  and  5.0  to  20.0  km/s  for  F+.  The  values 
found  in  the  present  study  agree  with  the  earlier 
findings,  albeit  on  the  low  side  of  the  range. 

V.  CONCLUSION 

The  species  present  in  a  PTFE  PPT  plume 
have  been  identified;  further,  the  composition  of 
the  PPT  exhaust  plasma  has  been  shown  to  be  in¬ 
dependent  of  discharge  energy  in  the  regions  sur¬ 
veyed.  Temporal  resolution  of  the  emission  has 
qualitatively  illustrated  the  dependence  of  emis¬ 
sion  on  discharge  energy  and  current.  Local  ther¬ 
modynamic.,  equilibrium  was  assumed  to  yield  an 
lectron  temperature  of  1.4+0. 2  eV  (~  16000  K). 
Spatially  and  temporally  resolved  images  were  ob¬ 
tained  to  place  bounds  on  the  axial  plasma  stream¬ 
ing  speed.  It  was  found  that  the  terminal  speed  of 
the  ionized  species  of  the  plasma  is  between  14  to 
16  km/s  after  the  first  half-cycle  of  the  discharge. 

More  accurate  quantitative  data  will  require 
spatial  resolution  of  the  discharge  plasma.  XPPT- 
1,  a  rectangular  geometry  thruster,  is  not  a  fa¬ 
vorable  design  on  which  to  perform  such  experi¬ 


ments.  However,  a  new  generation  of  coaxial  ge¬ 
ometry  thrusters  are  under  development  at  the 
USAF  Electric  _  Propulsion  Laboratory.  This  ax- 
isymetric  geometry  will  be  amenable  to  Abel  in¬ 
version  techniques  to  fully  spatially  resolve  the 
discharge  plasma  behavior. 


References 

[1]  J.  LeDuc,  D.  Bromaghim,  T.  Peterson, 
Mission  Planning.  Hardware  Development, 
and  Ground  Testing  for  the  Pulsed  Plasma 
Thruster  (PPT)  Space  Demonstration  on 
MightySat  II.  1  .  AIAA-97-2777,  33rd  Joint 
Propulsion  Conference,  Seattle  Wa.  July  7-9, 
1997. 

[2]  Cassady,  R.J.,  Myers,  R.M.,  and  Osborne, 
R.D.  Pulsed  Plasma  Thrusters  for  Spacecraft 
Attitude  Control.  JANNAF,  Albequerque, 
N.M.,  December,  1996. 

[3]  G.  Spanjers,  K.  McFall,  F.  Gulczinski  III. 
R.  Spores,  Investigation  of  Propellant  Ineffi¬ 
ciencies  in  a  Pulsed  Plasma  Thruster.  AIAA-' 
96-2723,  32nd  Joint  Propulsion  Conference, 
Lake  Buena  Vista,  FI.  July  1-3,  1996. 

[4]  Oriel  Corporation(1994),  Light  Sources, 
Monochromators,.  Spectrographs,  Detectors, 
Fiber  Optics  Vol.  II. 

[5]  Weast,  R.C.  CRC  Handbook  of  Chemistry 
and  Physics.  CRC  Press,  1989. 

[6]  Moore,  C.E.  Tables  of  Spectra  of  Hydrogen, 
Carbon,  ■  Nitrogen,  and  Oxygen  Atoms  and 
Ions.  CRC  Press,  1993. 

[7]  Pearse,  R.W.B.,  and  Gaydon.  A.G.  The 
Identification  of  molecular  Spectra.  Chapman 
and  Hall.  1976. 

[8]  Ruchti,  C.B.,  and  Niemeyer,  L.  Ablation 
Controlled  Arcs.  Transactions  on  Plasma  Sci¬ 
ence,  Vol.  PS-14,  No. 4,  August  1986. 

[9]  Auciello,  O.,  and  Flamm,  D.L.  Plasma  Di¬ 

agnostics  Volume  1.  Academic  Press.  Inc.. 
1989.  ‘  ' 

[10]  Demtroder,  W.,  Laser  Spectroscopy. 
Springer- Verlag,  1988. 

[11]  Bekefi.  G.,  Principals  of  Laser  Plasm, as.  Wi¬ 
ley,  1976. 

[12]  Vondra,  R..J.,  Thomasson,  K.I..  Solbes,  A., 
Exhaust  Velocity  Studies  of  a  Solid  Teflon 
Pulsed  Plasma  Thruster.  J.  Spacecraft  9,61 
(1972). 


10 


